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Abstract. The first science run of tlic LIGO and GEO gravitational wave detectors 
presented the opportunity to test methods of searching for gravitational waves from 
known pulsars. Here we present new direct upper limits on the strength of waves 
from the pulsar PSR J1939+2134 using two independent analysis methods, one in 
the frequency domain using frequentist statistics and one in the time domain using 
Bayesian inference. Both methods show that the strain amplitude at Earth from this 
pulsar is less than a few times 10~^^. 

PACS numbers: 04.80.Nn, 95.55.Ym, 97.60.Gb, 07.05.Kf 

Submitted to: Class. Quantum Grav. 

1. Introduction 

The LIGO and GEO gravitational wave detectors performed their first science run, 
denoted SI, from 23 August to 9 September 2002 (Abbott et al., 2003a, for the LIGO 
Scientific Collaboration). Although these instruments were still to reach their design 
sensitivities, their performances during this period were sufficiently good to justify 
a serious test of our search algorithms on real interferometer data. Code to search 
for continuous gravitational waves from rapidly rotating neutron stars has been under 
development within the LIGO Scientific Collaboration (LSC) since the mid-to-late 1990s, 
and for SI the LSC Pulsar Upper Limits Group (PULG) used two independent search 
methods to set upper limits on signals from the miUisecond pulsar PSR J1939-I-2134. 
This object has the shortest spin period of any known pulsar and falls in a relatively 
quiet part of the detector bands. Its short period also enables us to place a relatively 
tight upper limit on the equatorial ellipticity of the neutron star. Here we present a brief 
outline of the results from this exercise. A more detailed presentation of the methods 
and discussion of the results can be found in Abbott et al., 2003b (for the LIGO Scientific 
Collaboration) . 

2. Detector sensitivities 

In setting an upper limit on the gravitational wave amplitude from a pulsar we are 
saying something about the pulsar's physical parameters as constrained by instrumental 
performance rather than simply making a statement about the strain sensitivity of 
our instrument. The strength of any continuous gravitational wave signal from PSR 
J1939+2134 seen by a detecting interferometer will depend not only on the pulsar's 
distance, quadrupole moment and spin period, but also on the relative orientation 
between its axis of spin and the arms of the interferometer. The diurnal motion of 
the pulsar through the antenna pattern will also modulate the strength of the signal, so 
there is no simple relationship between the strain sensitivity curve of an interferometer 
and the upper limit that interferometer can place on the parameters of a particular 
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Figure 1. Mean upper limits, as a function of signal frequency, for the characteristic 
amplitude of gravitational waves from rotating neutron stars as constrained by the 
SI data (upper curves). The lower three curves represent ly design sensitivities for 
GEO, LIGO and one possible design for Advanced LIGO. The dotted lines show strain 
amplitudes corresponding to neutron stars at 8.5 kpc with equatorial ellipticities of 
10^'^, 10^^ and 10^^. The dots show characteristic amplitudes from known pulsars 
assuming all the observed loss in rotational kinetic energy is dissipated in gravitational 
waves. The arrow points to PSR J1939+2134. 



pulsar. We can however define a mean performance, averaged over tliese unknowns, for 
a particular data set and detection method. Figure [T] shows how our interferometers 
constrain monochromatic strain amplitudes at Earth as a function of signal frequency, 
averaged over the antenna pattern and relative orientation of the pulsar spin axis to the 
interferometer. The curves are derived from the 'frequency domain method' described 
below, and represent detection with a 1% false alarm rate and a 10% false dismissal rate 
using the actual observation times and data from SI. The lower solid curves show the 
design sensitivities of the interferometers for a 1 y observation period. The dotted lines 
indicate the strain amplitudes corresponding to hypothetical neutron stars at 8.5 kpc 
with equatorial ellipticities of 10~^, 10""^ and 10~^. Finally the dots represent the upper 
limit on strain amplitude for a selection of pulsars. For this we have assumed a simple 
spindown model in which all the apparent loss in rotational kinetic energy of the neutron 
star goes into the emission of gravitational waves. It is clear from this figure that no 
detection was expected from the SI data. Our primary aims in this exercise were to 
test the data reduction and analysis methods in the presence of realistic interferometer 
noise and set an upper limit on the strengths of such waves. 
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3. Seeirch methods 

Although the location, period and period derivative of PSR J1939+2134 are known to 
high precision the inclination of the neutron star's spin axis to the line of sight (i) and 
the phase (0), polarization (■0) and, of course, amplitude (ho) of the gravitational wave 
are unknown. All these factors affect the form of the signal and were searched over in 
SI using two independent methods. 

The 'frequency domain method' uses Fourier transform techniques to set a 
frequentist upper limit on the strain amplitude from the pulsar. The input time series 
is broken down into 60 s-duration short Fourier transforms (SFTs) which are high pass 
filtered at 100 Hz, Tukey windowed in the time domain and amplitude calibrated once 
per minute. These SFTs are then optimally combined to take account of the Doppler 
evolution and spindown of the signal using a Dirichlet kernel method (Jaranowski et al., 
1998) and used to compute a detection statistic, denoted J^, which is the likehhood ratio 
of the data, maximized over l, (f) and ip. Monte Carlo injections into the data stream 
are then used to calculate the probability distribution function (pdf) of J^. The upper 
limit is defined such that had a signal at the upper limit level or greater been present 
in an ensemble of similar experiments to the one performed, 95% of the experiments 
would yield values of that exceed the one seen in the real experiment. The computing 
requirement for this Monte Carlo work is extensive, and the searches ran off-line on the 
Medusa cluster at UWM, (296 single 1 GHz CPU nodes and 58TB of disk space), and 
the Merlin cluster at AEI Potsdam (180 dual 1.4 GHz CPU nodes and 36 TB of disk 
space) . 

The 'time domain method' is an algorithm specifically developed to search for 
neutron stars with a complex, but known, rotational phase evolution. It uses heterodyne 
methods to reduce and filter the data to a rate of one complex value every 60 s, estimating 
the noise variance over the same period to take account of non-stationarity in the 
data. The algorithm proceeds to carry out the search as a standard Baycsian modelling 
problem, determining a posterior pdf for the unknown source parameters ho, t, and 
i/j. The prior for is chosen to be uniform over [0, oo] and for ho < 0, and the three 
other priors are chosen to be least informative about the orientation of the neutron star. 
The time domain Bayesian upper limit is set via the posterior pdf for ho, marginalized 
over the nuisance parameters l, (j) and ip- We define the 95% upper credible limit for ho 
as the value at which the cumulative probability of ho (from zero) is 0.95. 

Note that there is nothing special about the choice of domain (frequency or time) 
for the frequentist and Bayesian analyses - the above choice was influenced by the desire 
to have two very different approaches to solving the same problem. 

4. Results 

As expected the SI data provided no evidence for continuous wave emission from PSR 
J1939-I-2134 at twice its rotation frequency (as determined by radio observations). 
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Instrument Frequentist UL Bayesian UL 



GEO 


(1.9 ±0.1) X 10- 


-21 


(2.2 ±0.1) 


X 10- 


-21 


LLO 


(2.7 ±0.3) X 10- 


-22 


(1.4 ±0.1) 


X 10- 


-22 


LH0-2k 


(5.4 ±0.6) X 10- 


-22 


(3.3 ±0.3) 


X 10- 


-22 


LH0-4k 


(4.0 ±0.5) X 10- 


-22 


(2.4 ±0.2) 


X 10- 


-22 



Table 1. Summary of the 95% upper limits delivered by the two analysis methods 
when applied to the four interferometers taking part in SI. GEO is the UK/German 
GEO 600 interferometer in Hanover, LLO is the LIGO instrument in Livingston, LA. 
LH0-2k and LH0-4k are the two LIGO interferometers co-sited at Hanford, WA. 



However, we were able to improve the best upper limit for such signals by about a 
factor of 100, based on earher work by Hough et al. (1983) and Hereld (1983). 

The resuhs from the four interferometers that took part in SI are summarized in 
Table Q the quoted uncertainties being largely due to calibration. The two analysis 
methods are very different, and the meaning of 'upper limit' is not the same for both, 
however we would hope that the concept of an upper limit is sufficiently robust that 
the two interpretations would deliver numerically similar results, as indeed they do. 
In principle the Bayesian analysis allows the data from all four interferometers to be 
combined to give an improved joint upper limit, although the relative phase calibration 
of the instruments was not sufficiently reliable to do this for SI. It is however interesting 
to note that an upper limit for ho can, within a suitable model, be interpreted as an 
upper limit on the equatorial ellipticity of the neutron star. If we take the distance to 
the pulsar to be 3.6 kpc, its radius to be 10 km and its mass to be 1.4 Mq a value of 
/lo < 1.4 X 10^^^ corresponds to an ellipticity of less than 2.7 x 10^^. Such a large upper 
limit to the ellipticity could be rejected a priori, as it implies a much greater spindown 
rate (via gravitational luminosity) than is seen in this pulsar. However as detector 
sensitivities improve we can expect the above methodology to deliver constraints on 
neutron star parameters not otherwise achievable. 
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